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ABSTRACT 

We present new, deep 8.5-GHz VLA observations of the nearby, low-luminosity ra- 
dio galaxy 3C296 at resolutions from 0.25 to 5.5arcsec. These show the intensity 
and polarization structures of the twin radio jets in detail. We derive the spectral- 
index distribution using lower-frequency VLA observations and show that the flatter- 
spectrum jets are surrounded by a sheath of steeper-spectrum diffuse emission. We 
also show images of Faraday rotation measure and depolarization and derive the ap- 
parent magnetic-field structure. We apply our intrinsically symmetrical, decelerating 
relativistic jet model to the new observations. An optimized model accurately fits the 
data in both total intensity and linear polarization. We infer that the jets are inclined 
by 58° to the line of sight. Their outer isophotes flare to a half-opening angle of 26° 
and then recoUimate to form a conical flow beyond 16 kpc from the nucleus. On-axis, 
they decelerate from a (poorly-constrained) initial velocity (3 = v/c~ 0.8 to /? « 0.4 
around 5 kpc from the nucleus, the velocity thereafter remaining constant. The speed 
at the edge of the jet is low everywhere. The longitudinal profile of proper emissiv- 
ity has three principal power-law sections: an inner region (0 - 1.8 kpc), where the 
jets are faint, a bright region (1.8 - 8.9 kpc) and an outer region with a flatter slope. 
The emission is centre-brightened. Our observations rule out a globally-ordered, helical 
magnetic-field configuration. Instead, we model the field as random on small scales but 
anisotropic, with toroidal and longitudinal components only. The ratio of longitudinal 
to toroidal field falls with distance along the jet, qualitatively but not quantitatively 
as expected from flux freezing, so that the field is predominantly toroidal far from 
the nucleus. The toroidal component is relatively stronger at the edges of the jet. A 
simple adiabatic model fits the emissivity evolution only in the outer region after the 
jets have decelerated and recollimated; closer to the nucleus, it predicts far too steep 
an emissivity decline with distance. We also interpret the morphological differences 
between brightness enhancements ("arcs") in the main and counter-jets as an effect of 
relativistic aberration. 

Key words: galaxies: jets - radio continuum: galaxies - magnetic fields - polarization 
- MHD 



1 INTRODUCTION 

The case that the jets in low- luminosity, FRI radio galaxies 
iFanarofF fc Rilevlll974) are initially relativistic and decel- 
erate to sub-relativistic speeds on kiloparsec scales rests on 
several independent lines of evidence, as follows. 

(i) Proper motions corresponding to speeds compara- 
ble with or exceeding c have been observed directly on 
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pc and kpc scales in FR I lets iGiovannini et al.l l200lt 
iBiretta. Zhou fc Owe Jl995l: iHardcastle et alJl2003^ . 

(ii) FRI sources must be the side-on counterparts of at 
least a subset of the BL Lac population, in which evidence 
for bu lk relativistic flow on pc scales is well established (e.g. 
lUrrv fc PadovanL.1995,,'1 . 

(iii) FRI jets show side-to-side asymmetries which de- 
crease with distance from the nucleus, most naturally ex- 
plained by Doppler beaming of emission from an intrinsi- 
cally symmetrical, decelerating, relativistic flow iLaing et alJ 
Il999h . in which case the brighter jet is approaching us. There 
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is continuity of sidedness from pc to kpc scales and superlu- 
minal motion, if observed, is in the brighter jet. 

(iv) The lobe containing the brighter jet is less depolar- 
ized than the counter-jet lobe, consistent with Faraday ro- 
tation in the surro unding halo of hot plasma if the b righter 
jet is approaching llLainelll988l : iMorganti et alJll997^ . 

(v) Theoretical work has demonstrated that initially 
relativistic jets can decelerate from relativistic to sub- 
relativistic speeds without disruption provided that they 
are not too powerful, the surrounding halo of hot plasma 
providing a pressure gradient large enough to rccoU imate 
them (Phinncv 1983; Bicknell . 1994; Komissarol Il994t 
iBowman. Leahv fc Komissarovlll99(i^l . 

We have developed a method of modelling jets on the as- 
sumption that they are intrinsically symmetrical, axisym- 
metric, relativistic, decelerating flows. By fitting to deep ra- 
dio images in total intensity and linear polarization, we can 
determine the three-dimensional variations of velocity, emis- 
sivity and magnetic-field ordering. We first a pplied the tech- 
nique to the jets in the radio galaxy 3C 31 iLaing fc Bridld 
l2002al. LB). With minor revisions, the same method proved 
capable of modelling the jets in B2 0326-1-39 and B2 1553-1-24 
(Canvin & Laine 2004, CL) and in the giant radio galaxy 
NGC315 dCanvin et ani2005l . CLBC). We used the derived 
velocity field for 3C31, together with Chandra observations 
of the hot plasma surrounding the parent galaxy, to de- 
rive the run of pressure, density, Mach number and en- 
trainment rate along its jets via a conservation-law analysis 
jLaing fc Brid le 20021^: 

We also examined adiabatic models 
for 3C31 in detaU iLaing fc Bridlell2004^ . 

We now present a model for the jets in 3C 296, derived 
by fitting to new, sensitive VLA observations at 8.5 GHz. 
In Section]^ we introduce 3C 296 and summarize our VLA 
observations and their reduction. Images of the source at 
a variety of resolutions are presented in Section |21 where 
we also describe the distributions of spectral index, rotation 
measure and apparent magnetic field derived by combining 
our new images with lower-frequency data from the VLA 
archive. We briefiy recapitulate our modelling technique in 
Section 2] and compare the observed and model brightness 
and polarization distributions in Section |3 The derived ge- 
ometry, velocity, emissivity and field distributions are pre- 
sented in Section |n] We compare the result of model fitting 
for the five sources we have studied so far in Section |7| and 
summarize our results in Section |S| 

We adopt a concordance cosmology with Hubble con- 
stant, Ho = 70kms~^Mpc~\ Q.a = 0.7 and Q.m = 0.3. 
Spectral index, a, is defined in the sense S{v) oc and 
we use the notations P — {Q^ + 11^)^''^ for polarized intensity 
and p = P/I for the degree of linear polarization. We also 
define /3 = v/c, where v is the fiow velocity. F = (1-/3^)"^''^ 
is the bulk Lorentz factor and 9 is the angle between the jet 
axis and the line of sight. 

2 OBSERVATIONS AND IMAGES 
2.1 3C296 

For our modelling technique to work, both radio jets must 
be detectable at high signal-to-noise ratio in linear polariza- 
tion as well as total intensity, separable from any surround- 



ing lobe emission, straight and antiparallel. 3C 296 satisfies 
all of these criteria, although the lobe emission is relatively 
brighter than in the other sources we have studied. This is 
a potential source of error, as we shall discuss. 

3C 296 was first image d at 1.4 and 2.7 GHz by 
iBirkinshaw. Laing fc Peacock! il98ll) . VLA observations 
showing the large-scale structure at 1.5 GHz and th e inner 
jets at 8.5 GHz were prese nted by iLeahv fc PerlevI (ll99J) 
and lHardcastle et al.l (Il997f) respectively. 

The large-scale radio structure of 3C 296 has two ra- 
dio lobes with well-defined outer boundaries and weak dif- 
fuse "bridge" emission aro und much of the pa th of both 
jets, at least in projection jLeahv fc Perlevlll99lll . 3C 296 is 
therefore mo re typical of a "bridged twin-jet" FRI structure 
jLainelll99el) than the "tailed twin-jet" type e xemplified by 
3C 31 and other sources we have studied so far. lPaxma et alJ 
( 1996) point out that the "lobed, bridged" FRI population is 
more abundant in a complete sample of low- luminosity radio 
galaxies than the "plumed, tailed" population. It is therefore 
particularly interesting to compare the inferred properties of 
the jets in 3C 296 with those of the sources we have modelled 
previously. Of the two lobes, that containing the counter-jet 
depolarizes more rapidly with decreasing frequency be tween 
1.7 and 0.6 GHz jGarrington. Holmes Sz SaikiallT99^ . con- 
sistent with the idea that it is further awa y from us. 

T he parsec-scale structure, imaged bv lGiovannini et alJ 
(|200^, is two sided, but brighter on the same side 
as the kpc-scale jet emission. Non-thermal X-ray 
and ultraviole t emission ha ve been detected from the 
brighter jet jHardcastle et a l. 2005) a nd the nucleus 
("Chiabcrg e. Capetti fc Celottil (1999t iHardcastle et all 
i2005). The source is associated with the giant elliptical 
galaxy NGC 5532 {z = 0.02470 ± 0.00007), the ce ntral, 
dominant member of a small group llMiUer et al.ll2002ll . The 
linear scale is 0.498 kpc arcsec"^ for our adopted cosmology. 
The galaxy has a slightly warped nuclear dust ellipse in 
PA 163° ± 4° with an inclination of 73° ± 4° assuming 
intrinsic circularitv l|Martel et aL|ll99St Ide'Koff et alJl2000l: 
IVerdoes Kleiin fcdeZeeu^ ~ l200g[ parameters are taken 
from the last reference). X-ray emission from hot plasma 
associated with th e galaxy was detected by ROSAT wad im- 
aged bv Chandra llMller et alJlTggl IHardcastle fc Worr'alil 
I1999I: IHardcastle et alJl2005^ . 

2.2 Observations 

New, deep VLA observations of 3C 296 at 8.5 GHz are pre- 
sented here. We were interested primarily in the inner jets, 
and therefore used a single pointing centre at the position 
of the nucleus and the maximum bandwidth (50 MHz) in 
each of two adjacent frequency channels. In order to cover 
the full range of spatial scales in the jets at high resolution, 
we observed in all four VLA configurations. We combined 
our new datasets with t he shorter observations described by 
IHardcastle et all il997l) . which were taken with the same 
array configurations, centre frequency and bandwidth. 

In order to correct the observed E-vector position an- 
gles at 8.5 GHz for the effects of Faraday rotation and to 
estimate the spectrum of the source, we also reprocessed 
earlier observations at 4.9 GHz and 1.4 - 1.5 GHz (L-band). 
The 4.9-GHz data were taken in the A configuration. At L- 
band, we used two A-configuration observations (the later 
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Table 1. Journal of VLA observations, v and Ap are the centre 
frequencies and bandwidth, respectively, for the one or two fre- 
quency channels observed and t is the on-source integration time 
scaled to an array with all 27 antennas operat ional. References 
to pre vio us publications using t he same data: l lHardcastle et alj 
<199'i) . 2 lLeahv fc PerlevI <199lD . 
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Date 


V 






t 


Ref 


uration 




(GHz) 




(MHz) 


(min) 




A 


1995 Jul 24 


8.4351, i 


J.4851 


50 


123 


1 


B 


1995 Nov 27 


8.4351, I 


J.4851 


50 


110 


1 


C 


1994 Nov 11 


8.4351, I 


5.4851 


50 


47 


1 


D 


1995 Mar 06 


8.4351, I 


5.4851 


50 


30 


1 


A 


2004 Dec 12 


8.4351, i 


5.4851 


50 


250 




A 


2004 Dec 13 


8.4351, I 


5.4851 


50 


253 




B 


2005 Apr 21 


8.4351, 1 


5.4851 


50 


250 




B 


2005 Apr 23 


8.4351, I 


5.4851 


50 


264 




C 


2005 Jul 14 


8.4351, I 


5.4851 


50 


249 




D 


2005 Nov 27 


8.4351, I 


5.4851 


50 


98 




A 


1981 Feb 13 


4.8851 




50 


130 




A 


1981 Feb 13 


1.40675 




12.5 


125 




A 


1995 Jul 24 


1.3851, 1.4649 


50 


45 


1 


B 


1987 Dec 08 


1.4524, 1.5024 


25 


60 


2 


C 


1988 Mar 10 


1.4524, 1.5024 


25 


50 


2 


D 


1988 Oct 7 


1.4524, 1.5024 


25 


80 


2 



one also from iHardcastle et alJ IlOQ?*! and the combined B, 
C and D-configuration dataset from Leahy & Pcrlev. LlQQLl . 
A journal of observations is given in Table Q 



2.3 Data reduction 

2,. 3.1 8.5 GHz observations 

All of the 8.5-GHz datasets were calibrated using standard 
procedures in the aips package. 3C 286 was observed as a 
primary amplitude calibrator and to set the phase difference 
between right and left circular polarizations; J1415-I-133 
was used as a secondary phase and amplitude calibrator 
and to determine the instrumental polarization. After ini- 
tial calibration, the 1994-5 data were precessed to equinox 
J2000 and shifted to the phase centre of the new observa- 
tions. All three A-configuration datasets were then imaged 
to determine a best estimate for the position of the core 
(RA 14 16 52.951; Dec. 10 48 26.696; J2000) and then self- 
calibrated, starting with a phase-only solution for a point- 
source model at the core position and performing one further 
phase and one amplitude-and-phase iteration with CLEAN 
models. The core varied significantly between observations, 
so we added a point source at the position of the peak 
to the 1995 dataset to equalize the core flux densities. We 
then combined the three A-configuration datasets and per- 
formed one further iteration of phase self-calibration. The 
B-configuration datasets were initially phase self-calibrated 
using the A-configuration model. After further iterations of 
self-calibration and adjustment of the core flux density, they 
were concatenated with the A-configuration dataset. The 
process was continued to add the C- and D-configuration 
datasets, with relative weights chosen to maintain a resolu- 
tion of 0.25 arcsec for images from the final combination. 
The off-source noise levels for the final I images (~ 4/iJy 



beam~^ rms at full resolution) are close to the limits set 
by thermal noise, although there are slight residual arte- 
facts within a pair of triangles centred on the core and de- 
fined by position angles ±15°. These do not affect the jet 
emission, and the maximum increase in noise level is in any 
case only to ~ 6/iJy beam~^ rms at full resolution. I im- 
ages were made using clean and maximum-entropy decon- 
volution, the latter being convolved with the same Gaus- 
sian beam used to restore the CLEAN images. Differences 
between the two deconvolutions were very subtle. As usual, 
the maximum-entropy method gave smoother images, avoid- 
ing the mottled appearance produced by CLEAN in regions 
of uniform surface brightness. On the other hand, the CLEAN 
images had slightly lower off-source noise levels, significantly 
lower (and effectively negligible) zero-level offsets, and im- 
proved fidelity in regions with high brightness gradients. 
The differences between the two methods are sufficiently 
small that none of our quantitative conclusions is affected 
by choice of method. We have used the CLEAN deconvolu- 
tions for all of the modelling described in Sections |1] - |S1 
but have generally preferred the maximum entropy images 
for the grey-scale displays in Section |3| All Q and U images 
were CLEANed. 

We made /, Q and U images at standard resolutions 
of 5.5, 1.5, 0.75 and 0.25 arcsec FWHM. All images were 
corrected for the primary-beam response during or after de- 
convolution, so the noise levels (listed in Table |2J strictly 
apply only in the centre of the field. The total flux density 
determined by integrating over a CLEAN image at 5.5 arcsec 
resolution is 1.31 Jy (with an estimated error of 0.04 Jy due 
to calibration uncertainties alone), compared with an ex- 
pected value of 1.20 ± 0.06 Jy, interpolate d between single- 
dish measurements at 10.7 and 5.0 GHz iLaing fc Peacock 
Il980l) . At least in the low-resolution images, we have there- 
fore recovered the total flux density of the source, despite 
the marginal spatial sampling (Table 



2.3.2 4.9 and 1.4-1.5 GHz observations 

For the B, C and D-configuration observations at 1.45 and 
1.50 GHz, we started from the combined dataset described 



bv lLeahy fc PerlevI (Il991f) . We imaged a significant fraction 
of the primary beam using a faceting procedure in order to 
minimize errors due to non-coplanar baselines and confusion, 
using the resulting CLEAN model as the input for a further 
iteration of phase-only self-calibration. This led to a signif- 
icant improvement in o fl-source noise level co mpared with 
the values obtained by iLeahv fc PerlevI Il99ll using earlier 
algorithms. The noise level in I is much larger at 1.45 GHz 
than at 1.50 GHz, probably as a result of interference; the 
Q and U images in the two frequency channels have simi- 
lar noise levels, however. We used both single-frequency and 
combined images at 5.5-arcsec resolution; the former for po- 
larization and the latter for spectral index. The reason was 
that the variation of Faraday rotation across the structure 
turned out to be large enough that the E- vector position an- 
gle difference between the two frequency channels changed 
signiflcantly across the structure (Section 13.31 . whereas the 
variations of spectral index did not lead to detectable differ- 
ences between the I images in the two channels (Section l3.2ll . 
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Table 2. Resolutions and noise levels for the images used in this 
paper, aj is the ofF-source noise level on the I image; crp the av- 
erage of the noise levels for Q and U . M and C denote maximum- 
entropy and CLEAN / deconvolutions. All noise levels were deter- 
mined before correction for the primary beam response and larger 
values are appropriate at significant distances from the field cen- 
tre. The app roximate maximum scale of struc ture imaged reliably 
is also given <Tavlor. Ulvestad fc Perlevl2nn3) . Note that the 1.5- 
arcsec L-band image includes A-configuration data at 1.41 GHz 
and BCD configuration data at 1.45/1.50 GHz (see text). 
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4.8 
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40 
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35 


40 
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1.40675 


ABCD 


27 






900 


0.75 


8.4601 


ABCD 
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4.2 


4.1 
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0.40 


8.4601 


ABCD 




4.2 




180 


0.40 


4.8851 


A 
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10 


0.25 


8.4601 


ABCD 
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The A-configuration observations at 1.39, 1.41, 1.46 and 
4.9 GHz were calibrated, self-calibrated and imaged using 
standard methods. The 1.41-GHz A configuration dataset 
had better spatial coverage (in multiple separated snap- 
shots) but poorer sensitivity compared with the later short 
observations at 1.39/1.46 GHz. We were able to make a sat- 
isfactory multi-configuration / image at 1.5 arcsec resolution 
using the 1.41-GHz dataset together with the 1.45/1.50-GHz 
BCD configuration combination. No core variability was ap- 
parent between observations. The frequency difference is suf- 
ficiently small that there should be no significant effects 
on the spectral index image of the jets (see Section 13.211 . 
The CLEAN deconvolution of this combined ABCD configu- 
ration L-band image image had significant artefacts, but the 
maximum-entropy image at this resolution was adequate, at 
least in the area of the jets. In order to determine rotation 
measures, however, the individual IFs were kept separate 
(exactly as in the lower-resolution case) and we therefore 
made Q and U images from the individual A-configuration 
datasets. 

The A-configuration image at 4.9 GHz shows only the 
innermost regions of the jets. No significant polarization is 
detected and we use it only to estimate the spectral index 
of the jets at high resolution. 

Details of the images are again given in Table |5| 
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Figure 1. Grey-scales of total intensity. A square-root transfer 
function has been used to emphasise the low-brightness emission, 
(a) 1.48-GHz BCD CLEAN image at 5.5-arcsec FWHM resolu- 
tion. The grey-scale range, — 80mjy (beam area)~^, is marked 
by the labelled wedge, (b) L-band ABCD maximum entropy im- 
age at 1.5-arcsec FWHM resolution from a combination of 1.41 
and 1.48-GHz datasets. The grey-scale range is - 5 mjy beam~^ . 
The areas covered by the two panels are identical. 
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Figure 2. (a) Grey-scale of total intensity for a maximum-entropy image of the jets of 3C 296 at 8.5 GHz. The resolution is 1.5arcsec 
FWHM and the grey-scale range is - 1.5 mjy beam~^. (b) Grey-scale of a Sobel-filtered version of this image. Narrow features in the 
brightness distribution ("arcs") are marked with their types as defined in Section lij.il 
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Figure 3. Grey-scales of total intensity for tlie jets of 3C 296 
at 8.5 GHz, both from maximum-entropy images, corrected for 
attenuation by the primary beam, (a) Image at 0.75-arcsec reso- 
lution. The grey-scale range is - 0.75 mjy beam~^ and the box 
shows the modelled area (Sectionsgl-EJ. (b) Image at 0.25 -arcsec 
resolution with a grey-scale range of — 0.4 mJy beam~^. 



3 TOTAL INTENSITY, SPECTRUM, 

FARADAY ROTATION AND APPARENT 
MAGNETIC FIELD 

3.1 Total-intensity images 

The overall structure of the source at 5.5-arcsec FWHM is 
shown in Fig. ^a)^. This image emphasises that the lobes 
are well-defined, with sharp leading edges. With the possible 
exception of the innermost ±15 arcsec, where the extended 
emission is very faint, the jets appear superimposed on the 
lobes everywhere. It is possible that the jets propagate (al- 
most) entirely within the lobes, although the superposition 
may also be a projection effect. Fig. 0^b), which shows the 
same area at a resolution of 1.5 arcsec FWHM, emphasises 
that the boundaries between the jets and lobes are not pre- 
cisely defined at 1.4 GHz, although the edges of the jets are 
marked by sharp brightness gradients. 

In contrast, the lobe emission is almost invisible in an 
8.5-GHz image at the same resolution (Fig.|2K)- The trans- 
fer function of this grey-scale has been chosen to emphasise 
sub-structure ("arcs") in the brightness distributions. These 
features are emphasised further and labelled in Fig. Hlb), 
which shows a S obel-filtered version of the same image. The 
Sobel operator (|Prattlll99 l'l computes an approximation to 
I V(/) I and therefore highlights large brightness gradients. 
Arcs are observed in both jets, but appear to have two char- 
acteristic shapes: (i) concave towards the nucleus and ap- 
proximately semicircular and (ii) oblique to the jet direction 
on either side of the jet but without large brightness gradi- 
ents on-axis. They do not occupy the full width of the jets. 
All of the arcs in the main jet are of type (i); of those in the 
counter-jet, the inner three marked in Fig.|2Ib) are of type 
(ii) but the outer two are of type (i). In Sectiou FG.?! we inter- 
pret the structural differences in the main and counter-jet 
arcs as an effect of relativistic aberration. 

Fig. 13 shows the inner jets of 3C 296 in grey-scale rep- 
resentations again chosen to emphasise the fine-scale struc- 
ture. The i nitia l "fiaring", i.e., decoUimation followed by rec- 
oUimation llBridld 119821) . is clearest at 0.75-arcsec resolu- 
tion (Fig. is the ridge of emission along the axis 
of the main jet noted bv iHardcastle et alJ (^93); the in- 
nermost arcs are also visible. Our highest-resolution image 
(0.25 arcsec FWHM; Fig. ^) shows the faint base of the 
main jet. This brightens abruptly at a distance of 2.7 arc- 
sec from the nucleus; further out there is complex, non- 
axisymmetric knotty structure within a well-defined outer 
envelope. All of these features are typical of FRI jet bases 
(e.g. LB, CL, CLBC). There is only a hint of a knot in the 
counter-jet at 1.8 arcsec from the nucleus, opposite the faint 
part of the main jet, but integrating over corresponding re- 
gions between 0.5 and 2.7 arcsec in the main and counter- 
jets gives /j//cj ~ 2.7. The counter-jet brightens at essen- 
tially the same distance from the nucleus as the main jet, 
but its prominent, non-axisymmetric, edge-brightened knot 
is significantly wider than the corresponding emission on the 
main jet side. 



This is essentially the same image as in figs 13 and 30 of 
iLeahv fc PerlevI <199ll) 
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Figure 4. False-colour plots of spectral index and intensity gradients, (a) Spectral index between 1.48 and 8.5 GHz at a resolution of 
5.5arcsec FWHM. The spectral-index range is 0.4-1.2. (b) Spectral index between 1.4 and 8.5 GHz at a resolution of 1.5arcsec FWHM. 
The spectral-index range is 0.4 - 0.9. In panels (a) and (b), values are plotted only if the rms error is <0.2; blanked pixels are coloured 
grey, (c) Sobel-filtered L-band image at 5.5arcsec FWHM; (d) Sobel-filtered 8.5-GHz image at 1.5arcsec FWHM. 



3.2 Spectrum 

The D-configuration observations at 8.5 GHz nominally sam- 
ple a largest spatial scale of only 180 arcsec. Although the 
total angular extent of 3C 296 is 440 arcsec, the accurate cal- 
ibration of our dataset and the fact that we have recovered 
the total flux of the source imply that we have enough infor- 
mation to construct a low-resolution image of spectral index, 
a, between 8.5 and 1.48 GHz. This is shown in Fig. 31 a). We 
have also made a spectral-index image at 1.5 arcsec resolu- 
tion between 8.5 and 1.41 GHz; here the values in the diffuse 
emission are not reliable, and we show only the inner jets 
(Fig. Iljj) . Sobel-filtered / images at the two resolutions are 
shown for comparison in Figs^Ic) and (d). 

Fig. BI^) shows that the low-brightness regions of the 
lobes have significantly steeper spectra than the jets. Spec- 
tral indices in the jets at 5.5-arcsec resolution range from 
a = 0.53 close to the nucleus in both jets to a ~ 0.6 
in the main jet and a ~ 0.8 in the counter-jet. The 
steepest-spectrum diffuse emission visible in Fig. UTa) has 
a « 1.3. As in other FRI objects jKatz-Stone fc RudniclJ 



119971: iKatz-Stone et al.lll999D . the jets appear to be super- 
posed on steeper spectrum lobe emission, so there must be 
blending of spectral components. A steeper-spectrum com- 
ponent is present on both sides of the main and counter-jets, 
although it is more prominent to the SE. It is clearly visible 
and resolved from « 17 arcsec outwards along the main jet. 
The apparent narrowness of the steep-spectrum rim within 
~50 arcsec of the nucleus in the counter-jet, is misleading, 
however. Much of the lobe emission detected at 1.5 GHz in 
this region (Fig. ^) is too faint at 8.5 GHz to determine a 
reliable spectral index and the image in Fig. BI^) is there- 
fore blanked. Conversely, any superposed diffuse emission 
will have a small effect on the spectral index where the 
jet is bright. Only within roughly a beamwidth of the ob- 
served edge of the jet are the two contributions comparable, 
with a total brightness high enough for a to be determined 
accurately. The spectrum therefore appears steeper there 
(Fig. BK)- With higher sensitivity, more extended steep- 
spectrum emission should be detectable around the inner 
counter-jet. 
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Figure 5. Profiles of tlie spectral index along the jets of 3C 296 in 
PA — 35.°86. Dashed line: a measured between 1.48 and 8.5 GHz 
at a resolution of 5.5 arcsec FWHM, along the jet ridge line from 
Fig.UTa'). Full line: o between 1.41 and 8.5 GHz at a resolution of 
1.5 arcsec, averaged over ±1.5"on either side of the jet axis, from 
Fig.Hb). 

Flatter-spectrum emission can be traced beyond the 
sharp outer bends in both jets, suggesting that the flow in 
the outer part of the radio structure remains coherent even 
after it deflects through a large angle, at least in projec- 
tion.^ A comparison between the spectral-index image and 
a Sobel-flltered L-band image (Figs and c) shows that 
the flatter-spectrum regions are bounded by high brightness 
gradients; they can therefore be separated both morpholog- 
ically and spectrally from the surrounding steep-spectrum 
emission. In the NE lobe, the flatter-spectrum "extension" 
of the jet can be traced until it terminates at a region with 
an enhanced intensity (and intensity gradient) on the South 
edge of the lobe (Figs UK and c). 

The spectral structure therefore suggests that the jets 
deflect when they reach the ends of the lobes rather than dis- 
rupting completely. In FR I sources such as 3C 296, no strong 
shocks (hot-spots) are formed, implying that the generalised 
internal Mach number of the flow, M = VP/TsPs < 1, 
where /3s c is th e internal sound speed and Fs — (1 — Ps)~^''^ 
For an ultrarelativistic plasma. Pa = 3 , 
so strong shocks need not occur at the bends even if the flow 
remains mildly relativistic (e.g. we infer a tr ansonic flow de- 
celerat ing from /3 0.8 to /3 « 0.2 in 3C 31; iLaing fc Bridie 
l2002d) . In 3C296, our models imply that P « 0.4 on-axis 
at a projected distance of 40 arcsec from the nucleus, with 
a possible deceleration to /? < 0.1 at larger distances (Sec- 
tion I^J. 

It is plausible that the spectra in the lobe regions 

^ Note, however, that the very flat spectral indices a < 0.5 seen 
at both ends of the source may be unreliable because of the large 
correction for primary beam attenuation required at 8.5 GHz. 
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Figure 6. (a) Normalized polarization gradient p'(0)/p(0) from a 
fit to p at frequencies of 8.5, 1.45 and 1.50 GHz over the range —20 
to -1-5 m~^. The resolution is 5.5 arcsec FWHM. (b) RM image 
at 5.5-arcsec FWHM resolution, from a fit to position angles at 
the frequencies used in panel (a), (c) RM at 1.5 arcsec from a fit 
to position angles at 8.5, 1.46 and 1.39 GHz. The labelled wedges 
for panels (b) and (c) are in units of radm~^. The areas covered 
by the panels correspond to Figs^for (a) and (b) and Fig. Olfor 
(c). Blanked pixels are coloured grey. 
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Figure 7. Profiles of RM and depolarization along the jet axis, 
taken to be in PA —35°. 86. Positive and negative distances refer 
to the main and counter-jet sides, respectively. Panels (a) - (c) 
show averages for boxes of length 60arcsec (along the jet axis) 
and width 39arcsec at a resolution of 5.5arcsec FWHM. A ±7.5- 
arcsec region around the core is excluded, (a) Mean normalized 
polarization gradient (p'(0)/p(0)}. (b) Mean RM. The horizontal 
line shows the estimated Galactic RM. (c) Rms RM. The full 
lines show ctrm after a first-order correction for fitting error, the 
dotted lines show the uncorrected values o-RMraw a^nd the curve 
shows a simple model fitted by eye to the data (see text), (d) and 
(e) As (b) and (c), but at a resolution of 1.5arcsec FWHM. The 
boxes are 18 arcsec long and are extended transverse to the jet to 
include all of the unblanked emission in Fig. inic). The excluded 
region around the core is ±1.5 arcsec. 



of 3C 296 have been steepened by a combination of syn- 
chrotron, inverse Compton and adiabatic losses but it is not 
clear from the spectral-index image whether the jets and 
their extensions propagate within the steep-spectrum lobes 
(as in the standard model for FRII sources) or appear su- 
perimposed on them. 

At higher resolution (Fig. |1]d) , there are variations in a 
both along and transverse to the jets. Fig.|H|shows a longitu- 
dinal profile of spectral index along the axis within 60 arcsec 
of the nucleus, where the superposed diffuse, steep-spectrum 
emission has a negligible effect. There is good agreement be- 
tween the spectral indices measured at 1.5 and 5.5-arcsec 



resolution in this area. The spectral index close to the nu- 
cleus in both jets is «0.62. The spectrum then flattens with 
distance from the nucleus to a ~ 0.53 at 40 arcsec. Further 
from the nucleus, the spectrum steepens again (Fig. ^p), 
but lobe contamination becomes increasingly important and 
it is not clear that we are measuring the true spectrum of 
the jet. The flattening with increasing distance in the first 
40 arcsec of both jets is in the opposite sense to that ex- 
pected from contamination by lobe emission, however. The 
few other FRI jets observed with adequate resolution and 
sensitivity (none of which are superposed on lobes) show 
a similar spectral flattening with distance: PKS 1333—3 3 
jKilleen et alJll98d). 3C449 llKatz-Stone fc Rudnick|[l997l) . 
NGC315 llLaing et all20o3) . and NGC 326 (Murgia, private 
communication) . 

A comparison of Figs|ljb) and (d) shows that the spec- 
tral index is flatter in both jets roughly where the more 
distinct arcs are seen, but the arcs are not recognisable in- 
dividually on the spectral index images. They contain only 
a few tens of percent of the total jet emission along their 
lines of sight, so a small difference in spectral index be- 
tween them and their surroundings would not be detectable. 
The apparent narrowness of the steeper-spectrum (a ~ 0.9) 
rim observed at the edges of the jets in Fig. Bib) results 
from the superposition of jet and lobe emission, just as at 
lower resolution. This effect confuses any analysis of intrinsic 
transverse spectral gradien ts in the jets, such as those in the 
flaring region of NGC 315 llLaing et al.ll2006h . but there are 
hints that the intrinsic spectrum flattens slightly away from 
the axis within «30 arcsec of the nucleus in 3C 296 (Fig.^Ja). 

At 0.4-arcsec resolution, we measure a « 0.67 between 
4.9 and 8.5 GHz for the bright region of the main jet between 
2 and 8 arcsec from the nucleus and a ~ 0.57 for the bright 
knot at the base of the counter-jet. 

3.3 Faraday rotation and depolarization 

Comparison of lower-resolution images at 1.7 and 0.6 GHz 
showed that the SW (counter-jet) lobe depolarize s more 
rapid ly than the NE (main jet) lobe (Garrin gton et alJ 
[19965. We derived the polarization gradient p'(0) = 
dp/d{X^) and the polarization at zero wavelength, p{0) from 
a linear flt to the degree of polarization as a function of 
A^, p(A^) ~ p(0) +p'{0)X^, using 5.5-arcsec resolution im- 
ages at 8.5, 1.50 and 1.45 GHz. The linear approximation 
is adequate for the low depolarization seen in 3C 296 and 
allows us to use images at more than two frequencies to re- 
duce random errors. We show an image of p'{0) /p(0) (closely 
related to depolarization) in Fig. 1^1 a). There is a gradient 
along the jets, in the sense that the outer counter-jet is more 
depolarized than the main jet; we show a profile along the 
axis in Fig. |3a). Both lobes show significant depolariza- 
tion: {p'(0)/p(0)) = -3.7 m"^ and -5.8 m"^ for the NE and 
SW lobes, respectively, corresponding to depolarizations of 
0.84 and 0.75 at 1.45 GHz. As can be seen from Figs|SJa) 
and|7{a), the diffuse emission in the NE lobe is significantly 
more depolarized ((p'(0)/p(0)) — — 4.8m~^) than the main 
jet. In the SW lobe, the depolarization of the diffuse emis- 
sion ((p'(0)/p(0)) — — 6.2m~^) is comparable to that of the 
counter-jet measured further than 100 arcsec from the nu- 
cleus (Fig.[7K). 

The integrated Faraday rotation mea- 
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Figure 8. Vectors with lengths proportional to the degree of polarization at 8.5 GHz and directions along the apparent magnetic field, 
superimposed on a grey-scale of the total intensity gradient (Sobel-filtered image) at 1.5 GHz. The resolution is 5.5arcsec FWHM. The 
vector directions are derived from the 3-frequency RM fits where there is adequate signal-to-noise ratio; elsewhere, the 8.5-GHz position 
angles are corrected using the mean RM over the area of the image. Vectors are plotted where P 2(jp and / baj. 



sure (RM) for 3C 296 is " 3 ± 2radm-2 
jSimard-Normandin. Kronberg fc Button! Il98ll) and a 
naean value of Oradm"'^ with an rms disper si on o f 
4radm~^ was derived by iLeahv. Poolev fc RilevI il98(jl 
from three- frequency imaging at a resolution of 6 x 22 arcsec^ 
FWHM. We imaged the distribution of RM at 5.5 arcsec 
resolution, fitting to E-vector position angles, Xi 8.5, 
1.50 and 1.45 GHz. Given the closeness of the two lower 
frequencies, we have poor constraints on the linearity 
of the X ~ relation, but the values of RM are suffi- 
ciently small that we can be sure that there are no nn 
ambiguities (this is confirmed by comparison w i th th e 
2.7-GHz images of iBirkinshaw. Laing &: Peacock! Il98ll) . 
The resulting RM image is shown in Fig. |5Jb). The mean 
RM is — 0.4radm~^, with an rms of 6.1radm~^ and a 
total spread of ±20radm~'^. Profiles of mean and rms 
RM along the jet axis are shown in Figs Elb) and (c). 



We made a first-order correction to the rms RM, cTRMraw, 
by subtracting the fitting error am in quadrature to give 
o"RM ~ (""RMraw ~ '^fitY^'^ ■ Both raw and corrected values 
are plotted in Fig.|7[c). There is a clear large-scale gradient 
in RM across the counter-jet lobe, roughly aligned with the 
axis and with an amplitude of ~8radm~^. Fluctuations 
on a range of smaller scales are also evident. The rms RM 
peaks close to the nucleus and decreases by a factor of 1.5 
- 2 by 200 arcsec away from it (Fig. |7J;). 

At 1.5-arcsec resolution, we made 3-frequency RM fits 
to images at 1.39, 1.46 and 8.5 GHz (Fig.l^J:)^. Note that the 
largest scale sampled completely in this image is «40 arcsec. 
Profiles of mean and rms RM at this resolution are shown in 

A position-angle image at 1.41 GHz (derived using only the A- 
configuration data at this frequency) is consistent with this RM 
fit but too noisy to improve it 
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Figure 9. Vectors with lengtlis proportional to tlie degree of po- 
larization at 8.5 GHz and directions along the apparent magnetic 
field, superimposed on a grey-scale of total intensity at the same 
frequency. A square root transfer function has been used to em- 
phasise the diffuse emission around the main jet. The resolution 
is 5.5 arcsec FWHM. The vector directions are derived from 3- 
frequency RM fits where there is adequate signal-to-noise ratio; 
elsewhere, the 8.5-GHz position angles are corrected using the 
mean RM over the area of the image. Vectors are plotted where 
P ^ 2(7p and / ^ Scr/. 



Figs Eld) and (e). There are no systematic trends in either 
quantity within 100 arcsec of the nucleus. 

Our data are quahtat ively consistent with the results 
of iGarrington et alJ Jl996l) in the sense that we see some- 
what more depolarization in the counter-jet lobe at 5.5- 
arcsec resolution. The corresponding profile of RM fluctua- 
tions is quite symmetrical, however. If we c alculate the depo- 
larizat ion expected with the beam used bv lGarrinerton et al.l 
il996h between 1.7 and 0.6 GHz using the observed degree of 
polarization at 1.4 GHz and the RM image at 5.5-arcsec res- 
olution, we flnd no gross differences between the lobes. This 
analysis neglects spectral variations across the source and 
depolarization due to any mechanism other than foreground 
Faraday rotation on scales >5.5 arcsec, however. 

A simple model in which the variations in RM result 
from foreground field fluctuations in a spherically-symmetric 
hot gas halo is reasonable for 30 296: the model profile plot- 
ted in F ig. (T^c) was derived following the prescription in 
iLaing et a l. (2006), assuming an angle to the line of sight 
of 58° (Section 16.211 and fitted by eye to the data. For the 
derived core radius = 300 arcsec, the asymmetries in the 
predicted ctrm profile for angles to the line of sight > 50° are 
small and reasonably consistent with our observations. A full 



analysis must, however, include the gradient of RM across 
the counter-jet lobe and the enhanced depolarization seen 
in that region, neither of which follow the same profile as 
the rms RM. The large-scale gradient is extremely unlikely 
to be Galacti c, as the total Galacti c RM estimated from 
the models of iDineen fc Colej ^2005^ is only -2.4radm"^ 
and fluctuations on 200-arcsec scales are likely to be even 
smaller. The gradient must, therefore, be local to 3C 296 (a 
similar large-scale RM gradient, al beit of much larger am - 
plitude, was observed in Hydra A bv lTavlor fc Perlevll993) . 
We defer a quantitative analysis of the Faraday rotation and 
depolarization until XMM-Newton imaging of the hot gas on 
the scale of the lobes is available for 3C 296. 

3.4 Apparent magnetic field 

We show the direction of the apparent magnetic field, ro- 
tated from the zero-wavelength E-vector position angles by 
90°. Over much of the observed region, this can be derived 
directly from the X ~ ^ fit. The 8.5-GHz images generally 
have more points with significant polarization than those at 
lower frequencies, however. At points with significant polar- 
ized signal at 8.5 GHz but no RM measurement, we deter- 
mined the apparent magnetic-field direction using the ob- 
served 8.5-GHz position angle and the mean RM for the 
region. The maximum error introduced by this procedure is 
~ 1.4° provided that the RM distribution is no wider than 
in Fig.lHJb). The apparent field directions are shown for the 
whole source at 5.5 arcsec FWHM resolution in Fig. |H1 in 
more detail for a small region around the nucleus at the same 
resolution in Fig. |5] and for the jets at 1.5-arcsec FWHM 
resolution in Fig. 1101 The vector lengths in these plots are 
proportional to the degree of polarizat ion at 8.5 GHz, cor- 
rected to first order for Ricean bias dWardle fc Kronberd 
1I97J); this is equal to p(0) within the errors at all points 
where we have sufficient signal-to-noise to estimate depolar- 
ization (Section l3.3ll. 

As noted bv lLeahv fc PerlevI il99ll) , the apparent mag- 
netic field in the lobe emission is circumferential, with a 
high degree of polarization at the edges of the lobes. The 
magnetic field directions are well aligned with the ridges of 
the steepest intensity gradients, as shown by the superpo- 
sition on the Sobel-filtered L-band image (Fig. (SJ. The jets 
show mainly transverse apparent field at 5.5 arcsec resolu- 
tion, but the low-brightness regions around them do not all 
exhibit the same pattern of polarization. Those to the NW 
of both jets show p > 0.4 with the field ordered parallel to 
the jet axis. Those to the SB of the main jet are comparably 
polarized but with the magnetic field oblique to the jet axis 
while those to the SE of the counter-jet are only weakly po- 
larized. Both lobes contain polarized structure related to the 
spectral- index structure evident in Fig.^Ia). The field struc- 
ture in the S lobe is suggestive of a fiow that turns through 
180° but retains a transverse-field configuration after the 
bend. That in the N lobe is suggestive of a jet extension 
whose magnetic configuration changes from parallel imme- 
diately after the bend to transverse. The emission at the 
S edge of the lobe where there is a locally strong intensity 
gradient also exhibits a high [p > 0.6) linear polarization 
with the magnetic field tangential to the lobe boundary, i.e. 
transverse to the jet "extension" suggested by our spectral 
data. 
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At 1.5-arcsec resolution, the jet field structure is re- 
vealed in more detail fFig. IIOII . Longitudinal apparent field 
is visible at the edges of both jets, most clearly close to the 
nucleus and at ~60arcsec firom it. There are enhancements 
in polarization associated with the two outermost arcs in 
each jet: in all four cases the apparent field direction is per- 
pendicular to the maximum brightness gradient. 

The 8. 5- GHz Q and U images at 0.75 and 0.25 arcsec 
FWHM used for modelling have been corrected for Faraday 
rotation using the 1.5-arcsec RM images interpolated onto 
a finer grid. The corresponding apparent field directions are 
discussed in Section 



4 THE MODEL 
4.1 Assumptions 

Our principal assumptions are as as in LB, CL and CLBC: 

(i) FRI jets may be modelled as intrinsically symmet- 
rical, antiparallel, axisymmetric, stationary laminar flows. 
Real flows are, of course, much more complex, but our tech- 
nique will still work provided that the two jets are statisti- 
cally identical over our averaging volumes. 

(ii) The jets contain relativistic particles with an energy 
spectrum n(E)dE = noE~^'^°'^^^ dE and an isotropic pitch- 
angle distribution, emitting optically-thin synchrotron radi- 
ation with a frequency spectral index a. We adopt a spectral 
index of a = 0.60 derived from an average of the spectral- 
index image in Fig.|ljb) over the modelled region. The max- 
imum degree of polarization is then po = (3a + 3) /(3q-|-5) = 
0.705. 

(iii) The magne tic fiel d is t a ngled on small scales, 
but anisotrop ic jLaind Il98lt I Begelman et alj Il984i 
iLaing fc BridlelE00 2a HLaing. Canvin fc Bridle 200fll . 



4.2 Outline of the method 

It is well known that synchrotron radiation from intrinsically 
symmetrical, oppositely directed bulk-relativistic outfiows 
will appear one-sided as a result of Doppler beaming. The 
jet/counter-jet ratio, 7j/7cj = /3 cos e')/(l - /3cos 6»)]^+" 
for a constant-speed, one-dimensional flow with isotropic 
emission in the rest frame. The key to our approach is that 
aberration also acts differently on linearly-polarized radia- 
tion from the approaching and receding jets, so their ob- 
served polarization images represent two-dimensional pro- 
jections of the magnetic-field structure viewed from different 
directions 6'-^ and 6'^-^ in the rest frame of the flow: sinSj = 
sin0[r(l - PcosO)]-^ and sinS^ = sin6i[r(l + /3cos0)]-^ 
This is the key to breaking the degeneracy between (5 and 
6 and to estimating the physical parameters of the jets 
JCanvin fc Laine 2004). 

We have developed a parameterized description of the 
jet geometry, velocity field, emissivity and magnetic-field or- 
dering. We calculate the emission from a model jet in /, Q 
and U by numerical integration, taking full account of rela- 
tivistic aberration and anisotropic rest-frame emission, con- 
volve to the appropriate resolution and compare with the 
observed VLA images. We optimize the model parameters 
using the downhill simplex method with (summed over 
7, Q and U) as a measure of goodness of fit (note that the 
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Figure 11. The geometry of the jet, showing the flaring and 
outer regions and the quantities which define the shape of its 
outer surface. 



modelling is not affected by Ricean bias as we fit to Stokes 
/, Q and U directly). The method is described fully by LB, 
CL and CLBC. 

Our method is restricted to the vicinity of the nucleus, 
where we expect that relativistic effects dominate over in- 
trinsic or environmental asymmetries. The modelled region 
for 3C 296 (indicated by the box on Fig. |^) is limited in 
extent by slight bends in the jets at ~45 arcsec from the 
nucleus. Contamination by surrounding lobe emission also 
becomes significant at larger distances. 



4.3 Functional forms for geometry, velocity, 
magnetic field and emissivity 

The functional forms used for velocity, emissivity and field 
ordering are identical to those described by CLBC except for 
a very small change to the emissivity profile (Section 14. 3. 41 . 
We briefly summarize the concepts in this Section and give 
the functions in full for reference in the Appendix (Ta- 
ble ETJ. 



J^.3.1 Geometry 

The assumed form for the jet geometry fFig. lllH is that used 
by CL and CLBC. The jet is divided into a flaring region 
where it first expands and then recoUimates, and a conical 
(almost cylindrical) outer region. The geometry is defined 
by the angle to the line of sight, 6, the distance from the 
nucleus of the transition between flaring and outer regions, 
ro, the half-opening angle of the outer region, and the 
width of the jet at the transition, xq (all fltted parameters). 

We assume laminar flow along streamlines defined by an 
index s (0 ^ s ^ 1) where s = corresponds to the jet axis 
and s = 1 to its edge. In the outer region, the streamlines are 
straight and in the flaring region the distance of a streamline 
from the axis of the jet is modelled as a cubic polynomial 
in distance from the nucleus. Distance along a streamline 
is measured by the coordinate p, which is exactly equal to 
distance from the nucleus for the on-axis streamline. The 
streamline coordinate system is deflned in Appendix 1X1 
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Figure 10. Vectors with lengths proportional to the degree of polarization at 8.5 GHz and directions along the apparent magnetic field, 
superimposed on a grey-scale of Sobel-filtered total intensity at the same frequency. The resolution is 1.5 arcsec FWHM. The vector 
directions are derived from 3-frequency RM fits where there is adequate signal-to-noise ratio; elsewhere, the 8.5-GHz position angles are 
corrected using the mean RM of the image. Vectors are plotted where P ^ 3a p and / ^ 5a j. The "arcs" are those marked in Fig. |5Jb). 



4.3.2 Velocity 

The form of the velocity field is as used by CL and CLBC. 
The on-axis profile is divided into three parts: roughly con- 
stant, with a high velocity /3i close to the nucleus; a linear de- 
crease and a roughly constant but lower velocity /3o at large 
distances. The profile is defined by four free parameters: the 
distances of the two boundaries separating the regions and 
the characteristic inner and outer velocities. Off-axis, the 
velocity is calculated using the same expressions but with 
truncated Gaussian transverse profiles falling to fractional 
velocities varying between Vi and vq fTable lXTl . 



(toroidal, orthogonal to the streamline in an azimuthal di- 
rection). We initially included a radial field component in 
the model, exactly as for NGC 315 (CLBC). The best-fitting 
solution had a small radial component on the jet axis, but 
not elsewhere, and the fitted coefficients were everywhere 
consistent with zero. The improvement in resulting from 
inclusion of a radial component was also very small (see 
Section [4.411 . We therefore present only models in which the 
radial field component is zero everywhere. The longitudi- 
nal/toroidal field ratio varies linearly with streamline index 
as in CLBC (Table ICT . 



4.3.3 Magnetic field 

We define the rms components of the magnetic field to be 
(Bf)^^^ (longitudinal, parallel to a streamline) and (Bt)^^^ 



4.3.4 Emissivity 

We take the proper emissivity to be eh, where e is the emis- 
sivity in / for a magnetic field B = {Bf + Bf)^^'^ perpen- 
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Figure 12. Comparison between data and model at a resolution of 0.75arcsce. The displayed region is ±40.95 arcsec from the rmcleus 
along the jet axis and the seale is indieated at the bottom of the Figure, (a) Observed eontours. The levels are —1, 1, 2, 4, 8, 16, 32, 64, 
128, 256, 512, 1024, 2048, 4096 x 20/iJy beam^^. (b) Model contours, with the same levels as in panel (a), (c) Profiles of total intensity 
along the jet axis. Full line: data; dashed line: model, (d) Grey-scale of observed jet/counter-jet sidedness ratio, in the range 0-3, chosen 
to show the variations in the outer parts of the modelled region, (d) Grey-scale of model sidedness ratio, (e) Profiles of sidedness ratio 
along the jet axis. Pull line: data; dashed line: model. 



dicular to the line of sight and h depends on field geometry: 
for I, ^ h ^ 1 and for Q and U — po ^ ^ -l-po. We 
refer to e, loosely, as 'the emissivity'. For a given spectral 
index, it is a function only of the rms total magnetic field 
and the normalizing constant of the particle energy distri- 
bution, e oc noB^'^". 

As in CL and CLBC, the majority of the longitudinal 
emissivity profile is modelled by three power laws e oc p^^ , 
with indices E = Sin, Ej^id and Eout (in order of increas- 
ing distance from the nucleus). The inner emissivity region 
corresponds to the faint inner jets before what we will call 



the bnghtening point. The middle region describes the bright 
jet base and the outer region the remainder of the modelled 
area. We have adopted slightly different ad hoc functional 
forms for the short connecting sections between these regions 
in CL, CLBC and the present paper. For 3C 296, we allow 
a discontinuous increase in emissivity at the end of the in- 
nermost region, followed by a short section with power-law 
index Sknot, introduced to fit the first knots in the main 
and counter-jets. We adopt truncated Gaussian profiles for 
the variation of emissivity with streamline index. The full 
functional forms for the emissivity (with the names of some 
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variables changed from CLBC in the interest of clarity) are 
again given in Table \KT\ 

4.4 Fitting, model parameters and errors 

The noise level for the calculation of for I is estimated 
as l/y/2 times the rms of the difference between the image 
and a copy of itself reflected across the jet axis. For linear 
polarization, the same level is used for both Q and U. This 
is the mean of 1/ \/2 times the rms of the difference image 
for Q and the summed image for U, since the latter is anti- 
symmetric under reflection for an axisymmetric model. This 
prescription is identical to that used by LB, CL and CLBC. 
The region immediately around the core is excluded from 
the fit. 

We fit to the 0.25-arcsec FWLIM images within 
8.75 arcsec of the nucleus and to the 0.75-arcsec FWHM im- 
ages further out. In doing so, we implicitly assume that the 
jet structure visible at 8.5 GHz is not significantly contami- 
nated by superposed lobe emission. For total intensity, this 
is justified by the faintness of the lobe emission close to the 
jet even at much lower resolution (Fig. [SJ; we discuss the 
possibility of contamination of the polarized jet emission in 
Section lS.3l In our initial optimizations, the quantities defin- 
ing the shape of the jet projected on the plane of the sky 
were allowed to vary. These were then fixed and the remain- 
ing parameters were optimized, using only data from within 
the outer isophote of the model. 

We derive rough uncertainties, as in LB, CL and CLBC, 
by varying individual parameters until the increase in 
corresponds to the formal 99% confidence level for inde- 
pendent Gaussian errors. These estimates are crude (they 
neglect coupling between parameters), but in practice give 
a good representation of the range of qualitatively reason- 
able models. The number of independent points (2588 in 
each of 3 Stokes parameters) is sufficiently large that we are 
confident in the main features of the model. The reduced 
Xred ~ 0.72, suggesting that our noise model is oversim- 
plified. A model with a radial field component and a more 
general variation of field-component ratios with streamline 
index has Xrad = 0.69, but the coefficients describing the 
radial field are all consistent with zero; in any case, the re- 
maining model parameters are all identical to within the 
quoted errors with those given here. 



5 COMPARISON BETWEEN MODELS AND 
DATA 

5.1 Total intensity 

The observed and modelled total intensities and jet/counter- 
jet sidedness ratios are compared in Figs 1121 - 1151 Figs 1121 
and 1131 show images and longitudinal profiles at 0.75 and 
0.25 arcsec resolution; averaged transverse profiles of total 
intensity and sidedness ratio are plotted in Figs 1141 and ll5l 
The following total-intensity features are described well 
by the model: 

(i) The main jet is faint and narrow within «3 arcsec of 
the nucleus; such emission as can be seen from the counter- 
jet is consistent with a similar, but fainter brightness distri- 
bution (Fig.inj. 



(a) Total intensity: data 




(b) Total intensity: model 
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Figure 13. Comparison between data and model at a resolution 
of 0.25 arcsec. The displayed region is ±9.625 arcsec from the nu- 
cleus along the jet axis and the scale is indicated at the bottom 
of the Figure, (a) Observed contours. The levels are —1, 1, 2, 4, 8, 
16, 32, 64, 128, 256, 512, 1024, 2048, 4096 x 20AtJybeam-i. (b) 
Model contours, with the same levels as in panel (a), (c) Profiles 
of total intensity along the jet axis. Full line: data; dashed line: 
model. 

(ii) At 3 arcsec, there is a sudden increase in brightness 
on both sides of the nucleus. 

(iii) The bright region at the base of the main jet extends 
to a distance of ~13 arcsec from the nucleus; except for the 
initial bright knot, the counter-jet is much fainter relative 
to the main jet in this region (Fig |12l :'l. 

(iv) The jet /counter-jet sidedness ratio also shows a ridge 
of nearly constant amplitude from 10 arcsec to the end of the 
modelled region (Figs ll2( l - f). The ratio varies from 2.2 on- 
axis to « 1 at the edge of the jet in the outer part of the 
fiaring region and in the outer region fFigs H^h . b). 

(v) The transverse total-intensity profiles from 15 arcsec 
to the end of the modelled region are accurately reproduced 
(Fig. [nil. 

5.2 Linear polarization 

Grey-scales and longitudinal profiles of the degree of polar- 
ization for model and data at a resolution of 0.75 arcsec are 
compared in Fig. 1161 and grey-scales of p at 0.25-arcsec res- 
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Figure 14. Comparison between averaged transverse profiles of 
total intensity in the main and counter-jets at a resolution of 
0.75 arcsec. (a) main jet, (b) counter-jet in the flaring region. The 
profiles are generated by averaging along along radii between 15 
and 27 arcsec from the nucleus and plotting against angle from 
the jet axis, (c) main jet, (d) counter-jet in the outer region. The 
profiles are generated by averaging along the jet between 27 and 
40.95 arcsec from the nucleus and plotting against distance from 
the jet axis. In all panels, the data arc represented by full lines 
and the model by dashed lines. 



olution are shown in Fig. 1171 The model grey-scales in the 
former figure are plotted for / > 5(J/ , to allow direct compar- 
ison with the observations, and without blanking, to show 
the polarization of the faint emission predicted at the edges 
of the jets. The observed polarization varies quite rapidly 
with distance from the nucleus in the flaring region, so a 
transverse profile of p averaged over a large range of dis- 
tances is not useful. We have, instead, plotted proflles at 
three representative positions, each averaged over 1.5 arcsec 
in distance from the nucleus fFigs [T5k - f). Since the po- 
sition angles vary significantly across the jet, we calculated 
scalar averages of p. In the outer region, the transverse p 
profile is almost independent of distance, but the polariza- 
tion at the edges of the jets in the outer region is difficult 
to see on the 0.75-arcsec resolution blanked images profiles, 
since individual points have low signal-to-noise ratios. In or- 
der to improve the accuracy of the transverse profiles we 
need to average over Q and U (as is done by the model fit) 
rather than the scalar p. To do this, we changed the origin 
of position angle to be along the jet axis. An apparent field 
parallel or orthogonal to the jet axis then appears entirely in 
the Q Stokes parameter (and we verified that U was indeed 
small). We then integrated Q and / along the jet axis and 
divided to give the profiles oi Q/I shown in Figs ll8f g'l and 
(h). The sign convention is chosen so Q > for a transverse 
apparent field and Q < for a longitudinal one. Finally, we 
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Figure 15. Comparison between averaged transverse profiles of 
jet/counter-jet sidedness ratio for data and model at a resolution 
of 0.75 arcsec. Full line: data; dashed line: model, (a) Averages 
along radii from the nucleus in the fiaring region between 15 and 
27 arcsec. (b) Average along the jet axis in the outer region be- 
tween 27 and 40.95 arcsec. 



plot vectors whose lengths are proportional to p and whose 
directions are those of the apparent magnetic field in Figs ll9l 
(0.75 arcsec) and 1201 (0.25 arcsec) with two different blank- 
ing levels as in Fig. 1161 All images of the observed degree 
of polarization have first-order corrections for Ricean bias 
iWardle & Kronbore 19741. 

The following features of the polarized brightness dis- 
tribution are well fitted by the model. 

(i) Both jets show the characteristic pattern of transverse 
apparent field on- axis with parallel field at the edges in the 
flaring region. 

(ii) The parallel-fleld edge of the main jet is visible 
throughout the bright region at high resolution, albeit with 
low signal-to-noise ratio (Fig. I20II . On-axis, the degree of 
polarization is significantly lower. 

(iii) At high resolution, the apparent field at the base of 
the counter-jet appears to be "wrapped around" the bright 
knot (Fig. . 

(iv) The longitudinal profile of p (Fig. 1161 1) is matched 
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Figure 16. Comparison between the degree of polarization, p = 
P/I, of tile data and model at a resolution of 0.75 arcsec. The 
displayed region is ±40.95 arcsec from the nucleus along the jet 
axis and the scale is indicated at the bottom of the Figure, (a) 
Grey-scale of observed degree of polarization, in the range - 0.7. 
Points are only plotted if the total intensity / > 5(t/ (Table |2j. 
(b) Grey-scale of the model degree of polarization in the same 
range, with the same blanking as in panel (a), (c) As in (b), but 
with no blanking applied, (d) Profiles of p along the jet axis. Full 
line: data (all points have / > 5it/, so no blanking has occurred); 
dashed line: model (also unblanked). 



well, reproducing the significant differences between the two 
jets. 

(v) Tfie transition between longitudinal and transverse 
apparent field (where p = 0) in tfie main jet is placed cor- 
rectly at a distance of 8 arcsec from the nucleus (Figs 1191 

im). 

(vi) In contrast, the counter-jet shows transverse appar- 
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(a) Degree of polarization: data 
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Figure 17. Comparison between the degree of polarization, 
p = P/I, of the data and model at a resolution of 0.25 arcsec. 
The displayed region is from —5 to +9.625 arcsec from the nu- 
cleus along the jet axis and the scale is indicated by the labelled 
bar. The cross shows the position of the core, (a) Grey-scale of 
observed degree of polarization, in the range - 0.7. Points are 
only plotted if the total intensity I > 5cr/ f Table 1^. (b) Grey- 
scale of the model degree of polarization with the same range and 
blanking as in panel (a). 



ent field on-axis even very close to the core (Figs 119b , and 
c). 

(vii) The asymptotic value of p at distances from the nu- 
cleus >30 arcsec is larger in the counter-jet (p ~ 0.5) than 
in the main jet {p ~ 0.4). 

(viii) The ridge of transverse apparent field is broader in 
the main jet than the counter-jet at distances from the nu- 
cleus >25 arcsec (Figs[TOli - c,[THk - h.lT^. 

(ix) The evolution of the transverse polarization profiles 
with distance from the nucleus in the flaring region, which 
differ considerably between the two jets, are modelled accu- 
rately (FigsElj, - f.lBl. 

(x) The apparent fleld at the edge of the outer region of 
the counter-jet is longitudinal, with a degree of polarization 
approaching the theoretical maximum of po ~ 0.7 (Fig.ll8h'). 



5.3 Critique 

Several features of the total-intensity distribution close to 
the nucleus are not well described by the model. The main 
problem is that the brightness distributions in both jets show 
erratic fluctuations causing the local jet/counter-jet sided- 
ness ratio to increase with distance from the nucleus. Our 
models postulate a decelerating flow, so the jet/counter-jet 
ratio is expected to decrease monotonically with distance. 
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Figure 18. Transverse profiles of the degree of polarization for 
the main jet (left) and counter-jet (right) at representative dis- 
tances from the nucleus. The resolution is 0.75arcsec FWHM. In 
all panels, the data are represented by full lines and the model 
by dashed lines, (a) - (f): flaring region. The profiles are gener- 
ated by taking scalar averages of p along radii from the nucleus 
and plotting against angle from the jet axis, (a) and (b): 15 - 
16.5arcsec from the nucleus; (c) and (d): 20.25 - 21.75 arcsec; (e) 
and (f) 25 — 26.5 arcsec. Panels (g) - (h): averages of Q/I along 
the jets between 27 and 40.95 arcsec from the nucleus, plotted 
against distance from the jet axis. 



The model for the main jet base also predicts slightly too 
narrow a brightness distribution (Fig. 1136 , and b): the fit to 
the longitudinal / profile of the main jet at high resolution 
(Fig. I13t ) is reasonable, but at lower resolution, where the 
jet is partially resolved, it appears worse (Fig. \12^']. As in 
all of the other sources we have modelled, 3C 296 shows de- 
viations from axisymmetry and erratic intensity variations 
in the bright region (Figsj^Jj and ll3ll : we cannot fit these. 

The counter-jet is slightly wider than the main jet 
in the flaring region: Fig. 1211 shows a comparison of the 
20/iJybeam~^ isophotes at 0.75 arcsec resolution. An in- 
evitable consequence is that /j//cj < 1 at the edges of the 
jet. This is shown most clearly in the averaged transverse 
sidedness profiles of Fig. 1151 In the flaring region, averaging 
between 15 and 27 arcsec from the nucleus, the minimum 
value of /j//cj ~ 0.8 on both sides. An intrinsically symmet- 
ric, outflowing, relativistic model cannot generate sidedness 
ratios <1, so there must be some intrinsic or environmental 
asymmetry. The outer region has an asymmetric transverse 
sidedness profile, with /j//cj < 1 on one edge only; this ap- 
pears to be caused by deviations from axisymmetry in the 
main jet (Figs 112b and l2H . 

Although the qualitative difference in width of the 
transverse apparent field region between the main and 
counter-jets is reproduced correctly, the model profile is still 
narrower than the observed one in the outer region of the 
main jet. The signal-to-noise ratio is low at the edges of 
the jets, and the model fits should not be taken too se- 
riously there, but the averaged profiles of Q/I shown in 
Fig. llSf g) suggest that there is a real problem: the longi- 
tudinal apparent field predicted by the model at the edges 
of the main jet (an inevitable consequence of the absence of 
a radial magnetic-field component) is not seen. In contrast, 
the corresponding profiles in the counter-jet are accurately 
reproduced (Fig. I18h ) . Between «20 and 60 arcsec from the 
nucleus, there is significantly more lobe emission at 8.5 GHz 
around the main jet than the counter-jet (Fig. ^p), so we 
have examined the possibility that superposition of unre- 
lated lobe emission causes the discrepancy between model 
and data. The diffuse emission around the main jet has an 
apparent field transverse to the jet axis, as required to ac- 
count for the discrepancy, but is not bright enough to pro- 
duce the observed effect. In order to test whether it is bright 
enough to affect the profile in Fig. I18f g'). we fit and sub- 
tracted linear baselines from the Q and / profiles for the 
jet and counter-jet, using regions between 15 and 20 arcsec 
from the axis as representative of the lobe emission. The 
subtracted profiles are indistinguishable from those plotted 
in Fig. I18f g) and (h) within 12 arcsec of the jet axis. Super- 
posed lobe emission cannot, therefore, be responsible for the 
discrepancy unless it is significantly enhanced in the vicinity 
of the jet. * 



In fact, the sidedness ratio has a low value at the position of 
the first knot in the counter-jet, 3 arcsec from the nucleus, 
compared with its typical value «10 between 4 and 9 arcsec: 
effectively, the first counter-jet knot is roughly a factor of two 
brighter than would be expected. The main jet base can be 
fit more accurately by a model with less low-velocity emis- 
sion close to the brightening point, but the counter-jet knot 
is then eliminated completely and the overall is increased. 



* There is significantly more confusing lobe emission at 1.5 GHz, 
as indicated by the steep-spectrum rim observed at 1.5-arcsec 
resolution fPig. HbV 
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(a) Data 




^b) Model 




(c) Un blanked model 




10 arcsec p = 1 



Figure 19. Vectors with lengths proportional to the degree of polarization, p, and directions along the apparent magnetic field, superim- 
posed on grey-scales of total intensity. The resolution is 0.75 arcsec and vectors are plotted every 0.6 arcsec. The polarization and angular 
scales are indicated by the labelled bars in the lowest panel and the areas plotted are the same as those in Figs ll2l and ll6l (a) data; (b) 
and (c) model. In panels (a) and (b), vectors are plotted only where I > 5(Tj (Table|5J; no blanking has been applied in panel (c). 
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Table 3. Fitted parameters and error estimates. 
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0.0 


0.0 
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Rough error estimates for were made by eye (see text). 



6 PHYSICAL PARAMETERS 

6.1 Summary of parameters 

In this section, all distances in linear units are measured 
in a plane containing the jet axis (i.e. not projected on the 
sky). The parameters of the best-fitting model and their 
approximate uncertainties are given in Table |21 

6.2 Geometry and angle to the line of sight 

The angle to the hne of sight from our model fits is 9 = 
58^2 de gree. The axis of the nuclear dust ellipse therefore 
differs from that of the jets b y ~ 15° in inclination and 
~ 31° in the plane of the sky iVerdoes Kleiin fc de Zeeuwl 



[200^. The orientation of the dust is such that the counter- 
jet would be on the receding side of the galaxy if the dust and 
jet axes are even ve ry roughly aligned. Note, however, that 
Iverdoes Kleiin fc de Zeeuw (^OO^) suggest that dust ellipses 
such as that seen in 3C 296 are preferentially aligned with 
the major axes of the galaxies and statistically unrelated to 
the jet orientations. 

The outermost isophote of the jet emission is well fitted 
by our assumed functional form. The jets flare to a maxi- 
mum half-opening angle of 26° at a distance of 8.3 kpc from 
the nucleus, thereafter recoUimating. The outer region of 
conical expansion starts at ro = 15.8 kpc from the nucleus 
and has a half-opening angle of ^0 ~ 5° . The outer envelope 
of the model jet emission is shown in Fig. I23f a) . Note that 
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(a) Data 




(b) Model 




Figure 20. Vectors with lengths proportional to the degree of 
polarization, p, and directions along the apparent magnetic field, 
superimposed on grey-scales of total intensity. The resolution is 
0.25 arcsec. The polarization and angular scales are indicated by 
the labelled bars in the lowest panel and the areas plotted are the 
same as those in Fig. 1171 (a) data; (b) and (c) model. In panels 
(a) and (b), vectors are plotted only where / > Scr/ (Table|3; no 
blanking has been applied in panel (c). 

our models fit to faint emission outside the 5a isophote of 
20/iJybeam~^ at 0.75-arcsec resolution. 

6.3 Velocity 

Contours of the derived velocity field are shown in Fig. 1221 
The positions of the velocity regions, with their uncertain- 
ties, and the profiles of velocity along the on-axis and edge 
streamlines are plotted in Figs 1231 b). (d) and (f). The ini- 
tial on-axis velocity is poorly constrained. The best-fitting 
value, /3i — 0.8, is very uncertain, and the upper limit is un- 
constrained by our analysis. The reason for this is that 
the fractional edge velocity vi = O.lto i is very low: the 
emission from the main jet close to the brightening point 




-O 



1 arcsec 



Figure 21. Contours of the 20/xJybeam~^ isophote for the main 
jet (heavy line) and the counter-jet (light line) rotated by 180° 
about the core. This figure illustrates that the counter-jet is 
slightly wider than the main jet in the flaring region. In the outer 
region, the two isophotes match well on one side, but not on the 
other (primarily because of a lack of symmetry in the main jet). 
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Figure 22. Contours of the model velocity field. The thin contour 
is at /3 = 0.005 and the thick contours are at intervals of 0.05 in 
the range /3 = 0.05 - 0.80. 



is dominated by slower material. Emission from any high- 
speed component is Doppler-suppressed and therefore makes 
very little contribution to the overall ■ In contrast, after 
the jets decelerate the on-axis velocity is well-constrained 
and consistent with a constant value of /3o = 0.40^qq3. 
The deceleration is extremely rapid in the best-fitting model 
(Fig. 1231 1) , but the uncertainties in the boundary positions 
are also large. The low fractional edge velocity at larger 
distances, vo = 0.04^q q2, is required to fit the transverse 
sidedness-ratio profile (Fig. 1151 and is well-determined. The 
fact that the sidedness ratio is less than unity in places im- 
plies that there is some intrinsic asymmetry (Section 15.311 . 
but the qualitative conclusion that the edge velocity is low 
remains valid. 

X-ray and ultraviolet emission is detected between 2 
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Figure 23. (a) The geometry of the jet, showing the flaring and outer regions, (b) and (c): sketches showing the relative positions of 
the boundary surfaces between velocity and emissivity regions. The boundaries are defined in Table IXTl and their positions for the best 
fitting model are given in Table ISl The full curves indicate the boundaries and the shaded areas their allowed ranges, also from Table ISl 
(b) Velocity. The regions of approximately uniform deceleration and asymptotic outer velocity are marked, (c) Emissivity. The region of 
enhanced emissivity between 1.8 and 9kpc is indicated, (d) and (e); profiles of intrinsic parameters along the jet axis in the rest frame, 
(d) the velocity profile, (e) no-B^^" derived from the emissivity, with no and B in SI units. Solid line: model; dashed line, adiabatic 
approximation with the magnetic-field structure expected from flux freezing. The latter curve is plotted twice, normalized to the model 
at 1.8 and 14kpc from the nucleus, respectively, (f) and (g): proflles of intrinsic parameters along the jet edge in the rest frame, (f) the 
velocity, (g) the emissivity encoded as in panel (e). 
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and ~10arcsec from the nucleus iHardcastle et al']l2005l) . 
corresponding to 1.2 - 5.9 kpc in the frame of the jet. This 
corresponds to the bright part of the flaring region, up to 
and possibly including the rapid deceleration. 

We have evaluated the jet/counter-jet sidedness ratio 
outside the modelled region by integrating flux densities 
in boxes of size lOarcsec along the jet axis and ±15arcsec 
transverse to it. The mean ratio drops from «1.7 between 
40 and 70arcsec from the nucleus to ~1.05 between 70 and 
100 arcsec. It is therefore possible that further deceleration 
to /3 < 0.1 occurs at ~40kpc in the jet frame. As noted 
in Section 14.21 however, the jets bend and lobe contamina- 
tion becomes significant, so we cannot model the velocity 
variation in this region with much confidence. 



6.4 Emissivity 

The positions of the emissivity regions, with their uncer- 
tainties, and the profiles of emissivity along the on-axis and 
edge streamlines are plotted in Figs l23f c'l. (e) and (g). As in 
all of the FR I sources we have examined so far, the on-axis 
emissivity profile of 3C 296 is modelled by three power-law 
sections separated by short transitions. Close to the nucleus, 
before the brightening point at 1.8 kpc, the power-law index, 
Ein ~ 2.5 ± 0.5. The errors on Ein are estimated by eye be- 
cause the jets are faint and any change to the index has very 
little effect on the value of x^, even over a small area close to 
the nucleus. We model the brightening point as a discontin- 
uous increase in emissivity by a factor of l/g ~ 7. The short 
region between 1.8 and 2.5 kpc is introduced primarily to 
model the first knots in the main and counter-jets; its index 
is poorly constrained because of our inability to fit the main 
and counter-jet simultaneously: high values fit the counter- 
jet knot well but overestimate the main jet brightness; small 
values have the opposite problem. The main section of the 
bright jet base from 2.5 - 8.9 kpc is described by a power 
law of index i?rnid = 2.8. This joins smoothly onto the outer 
region (-Bout = 1.0) without the need for a rapid transition. 
The transverse emissivity profile is fairly flat at the bright- 
ening point (albeit with large uncertainties) but evolves to 
a well-deflned Gaussian with a fractional edge emissivity of 
0.12 at large distances. 

6.5 Magnetic-field structure 

The inferred magnetic fleld components are shown as pro- 
flies along the on-axis and edge streamlines in Fig. 1241 and 
as false-colour plots in Fig. 1251 In the former figure, the 
shaded areas define the region which the proflle could oc- 
cupy if any one of the six free parameters deflning it is varied 
up to the error quoted in in Table El We were able to flt the 
polarization structure in 3C 296 accurately using a model 
with only toroidal and longitudinal magnetic-field compo- 
nents. The best fitting model has field-component ratios 
which vary from the nucleus out to a distance of 14.6 kpc, 
consistent with the boundary between fiaring and outer re- 
gions, and thereafter remain constant. The field at the edge 
of the jet is almost purely toroidal everywhere (Figs 124b . 
I25t i') . although there is a marginally significant longitudinal 
component close to the nucleus (^Figs l24h . I25b 'l. On-axis, 
in contrast, the longitudinal fleld is initially dominant, but 



decreases over the fiaring region to a smaller, but still sig- 
nificant value at larger distances (Figs l^lb . 125b '). 

We have described the magnetic field structure as disor- 
dered on small scales, but anisotropic. A field of this type in 
an axisymmetric model always generates symmetrical trans- 
verse profiles of total intensity, degree of polarization and 
apparent field position angle. In contrast, a globally-ordered 
helical field, unless observed at 90° to the line of sigh t in its 
rest frame, will always generate asy mmetric profiles ijLaind 
Il98ll: iLaing. Ca nvin & Bridlell2006h . The observed profiles 
in 3C 296 (Figs 1141 1181 are very symmetrical, even in re- 
gions where the longitudinal and transverse fleld components 
have comparable amplitudes, so we conclude that a simple, 
globally-ordered field configuration cannot fit the data. Our 
calculations would, however, be unchanged if one of the two 
field components were vector-ordered; in particular, a con- 
figuration in which the toroidal component is ordered but 
the longitudinal one is not would be entirely consistent with 
our results. 

The polarization in the high-signal regions closer to the 
jet axis requires a mixture of longitudinal and toroidal field 
components, implying in turn that p ~ po at the edge of the 
jet with an apparent field along the jet axis. Our models in- 
clude this emission, but at a level below the blanking criteria 
in total intensity and linear polarization at 0.75-arcsec reso- 
lution (Fig. 116b '). We demonstrated in Section f5.3l that this 
emission is indeed present and that it has the correct po- 
larization in the counter-jet. In the main jet, contamination 
by surrounding lobe emission confuses the issue, but there 
is a real discrepancy (Fig |18b ): the observed apparent field 
remains transverse at the edges. The difference between the 
main and counter-jets cannot be an effect of aberration as 
we infer velocities /? ^ 1 in these regions. There must, there- 
fore, be a significant radial field component at the edge of 
the main jet which is not present at the corresponding loca- 
tion in the counter-jet. Although the low-resolution 8.5-GHz 
image shown in Fig.|5]gives the impression that the jet might 
propagate within its associated lobe while the counter-jet is 
in direct contact with interstellar medium, the L-band im- 
ages (Fig. show diffuse emission surrounding both jets, so 
this is unlikely to be the case. 

6.6 Flux freezing and adiabatic models 

Given the assumption of fiux freezing in a jet without a 
transverse velocity gradient, the magnetic field components 
evolve according to: 

Bt oc (x/JF)"^ 
ni (X X 

in the quasi-one-d imensional approx imation, where x is the 
radius of the jet jBaum et all 1199711 . The field-component 
evolution predicted by these equations is shown by the 
dashed lines in Fig. 1241 arbitrarily normalised to the model 
values at a distance of 14 kpc from the nucleus. The evolu- 
tion of the longitudinal/toroidal ratio towards a toroidally- 
dominated configuration is qualitatively as expected for flux- 
freezing in an expanding and decelerating fiow, but the 
model shows a much slower decrease. In the absence of any 
radial field, the ratio should not be affected by shear in the 
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Figure 24. Profiles of magnetic-field components along the axis of the jet (panels a - b) and its edge (panels c - d). The solid lines show 
the best- fitting model, the shaded areas the uncertainties derived from the limits in Table l3l and the dashed lines the profiles expected 
for a magnetic field frozen into the fiow. The profiles for a passively convected field are normalized to the free model predictions at a 
distance of 14 kpc from the nucleus, (a) and (c) toroidal; (b) and (d) longitudinal. 



laminar velocity field we have assumed, but the field evo- 
lution is clearly more complex than such a simple picture 
would predict. A small radial field component could still 
be present, and this would allow the growth of longitudi- 
nal field via the strong shear that we infer in 3C 296. We 
will investigate this proce ss using the ax i symm etric adia- 
batic models developed bv lLaing fc BridT^ (l2004l) . A second 
possibility is that the fiow is turbulent in the fiaring region. 
Two-dimensional turbulence (with no velocity component in 
the radial direction) would lead to amplification of the longi- 
tudinal component without generating a large radial compo- 
nent, as required. After the jets recoUimate, the model field 
ratios have constant values on a given streamline, consis- 
tent with the variation expected from fiux freezing in a very 
slowly expanding fiow (Fig. I24t . Note, however, the prob- 
lems in fitting the main jet polarization at large distances 
from the axis (Section 16.511 . 

If the radiating electrons sufi'er only adiabatic losses, 
the emissivity is: 



using equation 8 of iLaing fc Bridi3 ||200' 



-{l + 2c,/3) j^l + a 



in the quasi-one-dimension al approximation jBaum et al.l 
Il997t iLaing fc Bridldliooi) . B can be expressed in terms 
of the parallel-field fraction / = {Bf)^^'^ / B and the radius 
X, velocity /3 and Lorentz factor F at some starting location 



B oc 



1/2 



The resulting profiles for the on-axis and edge streamlines 
are plotted as dashed lines in Fig. I23f e) and (g), normal- 
ized to the model values at the brightening point (1.8 kpc) 
and 14 kpc. Their slopes are grossly inconsistent with the 
model emissivity in the fiaring region, but match well be- 
yond ~ 14 kpc - essentially where the jet recoUimates at ro 
= 15.8 kpc. 



6.7 Arcs 

As pointed out in Section 13.11 there are discrete, narrow 
features ("arcs") in the brightness distributions of both jets 
in 3C 296. Strikingly similar structures are seen in the jets 
of 3C31 (LB), as shown in Fig. ^ In both sources, the 
arcs have systematically different shapes in the main and 
counter-jets: centre-brightened arcs with well-defined inten- 
sity gradients near the jet axis [type (i)] are found in the 
main jets of both sources and in the outer counter-jet of 
3C 296, while more elongated arcs with well-defined inten- 
sity gradients near the jet edges [type (ii)] are found in both 
counter-jets. We will present models of the arc structures 
elsewhere: here we demonstrate that relativistic aberration 
can plausibly account for the systematic difference. 
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Figure 25. False-colour images of the fractional magnetic- field 
components, (a) toroidal {B^y^'^ /B; (b) longitudinal {Bf)'^^^/B. 
B= {B^ + Bfy/'^. 



Suppose that an arc is a thin, axisymmetric shell of 
enhanced emissivity, concave towards the nucleus, and that 
at any point its speed is roughly that of the local average 
flow. Close to the centre-line of the jet, the shell can be 
approximated as a planar sheet of material orthogonal to 
the axis in the observed frame. Aberration causes a moving 
object to appear at the angl e to the line of sight approp riate 
to the rest frame of the flow llPenrosell95 a lTerrelllQSa) . For 
the the normals to the sheets in the main and counter-jets 
these angles are O'j and where: 

sm9j = sin6'[r(l -/Scos^)]"^ 
sin^^ = sine[r(l -f/Scos^)]"^ 

just as in Section f4. 21 

In 3C 296, the arcs are visible as far as ~90 arcsec from 
the nucleus - outside the region we model - so we do not have 
a direct measure of the average velocity field in their vicinity. 
Given that our model indicates a roughly constant velocity 
on-a:xis from 10 - 40 arcsec from the nucleus we will assume 
for the purposes of a rough estimate that the parameters at 
the outer edge of the model approximate the average fiow pa- 
rameters near the arcs between 40 and 70 arcsec, where the 



mean jet/counter-jet sidedness ratio remains significantly 
larger than unity. In Fig. 1271 we plot 9'j and d'^j against (3 
for e = 58° (Table EJ. For a velocity /3 = /3o = 0.40, the 
angles to the line of sight in the fiuid rest frame are 9'j = 80° 
and 9'^j = 40° (the main jet is close to the maximum boost 
condition /3 = cos 9 which corresponds to edge-on emission 
in the rest frame). A sheet of enhanced emissivity which 
is normal to the flow on the jet axis is therefore observed 
nearly edge-on in the main jet, appearing narrow, with a 
high brightness gradient. In the counter-jet, on the other 
hand, the same sheet would appear to be significantly ro- 
tated about a line perpendicular to the jet axis in the plane 
of the sky and therefore less prominent both in intensity 
and in brightness gradient. Close to the edges of the jets, 
however, the model flow velocities are essentially zero, so 
the effects of aberration should be negligible. If the shells 
are roughly tangential to the surface, then we expect to see 
narrow brightness enhancements orientated parallel to the 
edges in both jets. 

This picture provides a reasonable qualitative descrip- 
tion of the differences in arc structure for the main and 
counter-jets in 3C 296 within ~70 arcsec of the nucleus 
(Fig. |5J . The evolution of the mean jet /counter-jet ratio sug- 
gests that the jets become sub-relativistic at larger distances 
(Section 16.31 1 and this may explain why we see type (i) arcs 
in both jets between 70 and 90 arcsec. 

The hypothesis that the main difference in appearance 
of the arcs between the main and counter-jets is an effect 
of differential relativistic aberration also allows us to pre- 
dict that differences should be observed between main and 
counter-jet arcs in other FRI sources if they are orientated 
at moderate angles to the line of sight. This explains why we 
observe systematic differences in the arcs of 3C 31, which has 
similar model parameters to 3C 296 (Fig. 1261 LB), but we 
would not expect differences between the main and counter- 
jet arcs if 9'j and 9'^j are similar. The latter condition holds 
if the jets are close to the plane of the sky (e.g. 3C 449 and 
PKS 1333-33: iFeretti et aljIigigtlKilleen et al.ll986^ or slow 
(e.g. the outer region of B2 0326-1-39; CL). It will therefore be 
interesting to search for arcs in deep images of these sources. 

7 COMPARISON BETWEEN SOURCES 

We have modelled five FRI sources: 3C31 (LB), 
B20326+39, B2 1553-^24 (CL), NGC315 (CLBC) and 
3C296 (this paper). We will compare and contrast their 
properties in detail elsewhere, but in this section we sum- 
marize their similarities and differences and draw attention 
to a few potential problems with our approach. 

7.1 Similarities 

The five FRI sources we have studied in detail so far all 
exhibit the following similarities: 

(i) A symmetrical, axisymmetric, decelerating jet model 
fits the observed brightness and polarization distributions 
and asymmetries well. 

(ii) The outer boundaries of the well-resolved parts of the 
jets can be divided into two main regions: a flaring region 
with a rapidly increasing expansion rate followed by recoUi- 
mation and an outer region with a uniform expansion rate. 
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Figure 26. 8.5-GHz images of tlie jets in the FR radio galaxy 3C31 (LB) with grey-scale ranges chosen to emphasise the arc structures, 
as in the equivalent plots for 3C296 (Fig. 1^. (a) Total intensity, (b) Sobel-filtered /. The resolution is 0.75arcsec FWHM. The arcs are 
labelled with their types as defined in Section Ik. 71 (except for the fourth from the nucleus in the main jet, which appears to be associated 
with the bend in the jet and is clearly not axisymmetric). 



(iii) Close to the nucleus, the on-axis velocity is initially 
consistent with /3 « 0.8. On-axis, rapid deceleration then 
occurs over distances of 1 - 10 kpc, after which the velocity 
either stays constant or decreases less rapidly. 

(iv) All except B2 1553-1-24 have low emissivity at the 
base of the "flaring" region, so even the main jet is initially 
very faint (this region may be present but unresolved in 
B21553-f24). 

(v) The faint region ends at a "brightening point" where 
the emissivity proflle flattens suddenly or even increases. 
This, coupled with the rapid expansion of the jets, leads to 
a sudden increase in surface brightness in the main jets (not 
necessarily in the counter-jets, whose emission is Doppler 
suppressed) . 

(vi) The brightening point is always before the start of 



rapid deceleration, which in turn is completed before the 
jets recollimate. 

(vii) Immediately after the brightening point, all of 
the jets contain bright substructure with complex, non- 
axisymmetric features which we cannot model in detail. 

(viii) Optical/ultraviolet or X-ray emission has been de- 
tected in all of the jets observed so far with HST or Chandra 
(results for B2 0326-1-39 are not yet available). High-energy 
emission is associated with the bright regions of the main 
jets before the onset of rapid deceleration in all cases, and 
also with the faint inner region in 3C 31. 

(ix) The longitudinal emissivity profile may be repre- 
sented by three power-law segments separated by short tran- 
sitions. The indices of the power laws decrease with distance 
from the nucleus. 
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(x) All of the jets are intrinsically centre-brightened, with 
fractional emissivities at their edges «0.1 - 0.5 where these 
are well-determined. 

(xi) On average, the largest single magnetic field compo- 
nent is toroidal. The longitudinal component is significant 
close to the nucleus but decreases with distance. The radial 
component is always the smallest of the three. 

(xii) For those sources we have fit with functional forms 
allowing variation of field-component ratios across the jets 
(3C 31, NGC 315 and 3C 296), we find that the toroidal com- 
ponent is stronger relative to the longitudinal component at 
the edge of the jet (this is likely to be true for the other two 
sources) . 

(xiii) The evolution of the field component ratios along 
the jets before they recoUimate is not consistent with fiux 
freezing in a laminar flow, which requires a much more rapid 
transition from longitudinal to transverse field than we in- 
fer. A possible reason for this is two-dimensional turbulence 
(with no radial velocity component) in the flaring regions. 
After recollimation, flux freezing is consistent with our re- 
sults. 

(xiv) The inferred evolution of the radial/toroidal field 
ratio (in those sources where the radial component is non- 
zero) is qualitatively inconsistent with fiux freezing even if 
shear in a laminar velocity field is included. 

(xv) The quasi-one-dimensional adiabatic approximation 
(together with flux freezing) is grossly inconsistent with 
the emissivity evolution required by our models before the 
jets decelerate and recoUimate. Earlier attempts to infer 
fracti onal changes in velocity from surface-brightness pro- 
flles (Fanti et al. 'l982'; BickncU '1984, 'l986; "Bickn eU et aP 
fl99Q,; Baum et al. .1997; Bondi et al. 2000) , therefore re- 
quired much more rapid deceleration than we deduce in these 
regions. 

(xvi) After deceleration and recollimation, the emissiv- 
ity evolution along the jets is described quite well by fiux 
freezing and a quasi-one-dimensional adiabatic approxima- 
tion (our analysis does not cover this region in NGC 315). 
As relativistic effects are slight, the brightness evolution in 
the outer regions of 3C 296, B2 03264-39 and B2 1553-1-24 is 
quite close to the constant - speed, p erpe ndicular- field predic- 
tion originally derived bv lBurchI il979f) . Previous estimates 
of fractional velocity changes derived on the adiabatic as- 
sumption [references as in (xv), above] are likely to be more 
reliable on these larger scales. 



7.2 Differences 

3C 296 exhibits the following differences from the other 
sources: 

(i) The transverse velocity profile falls to a low fractional 
value < 0.1 at the edge of the jets in 3C 296, whereas all 
the other sources are consistent with an edge/on-axis veloc- 
ity ratio «0.7 everywhere, and also with a top-hat velocity 
profile at the brightening point. This could be related to the 
differences in lobe structure: as noted earlier, 3C 296 is the 
only source we have observed with a bridged twin-jet struc- 
ture. Its jets may then be embedded (almost) completely 
within the lobes rather than propagating in direct contact 
with the interstellar medium of the host galaxy. Modelling of 
other bridged FRI sources would then be expected to show 
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Figure 27. The angle to the line of sight in the fluid rest frame for 
a vector parallel to the flow in the observed frame plotted against 
j3 = v/c. The angle between the jet axis and the line of sight is 
9 = 58°, as we infer for 3C 296. Full line: main jet; dashed line: 
counter-jet. The vertical dotted line marks the on-axis velocity of 
P = 0.4 we derive for the outer parts of the modelled region in 
3C296. 



similarly low edge velocities. Whether the lower-velocity ma- 
terial is best regarded as part of the jets or the lobes (and, 
indeed, whether this distinction is meaningful) remains un- 
clear. 

(ii) Some support for the idea of an interaction between 
jets and lobes comes from our detection of anomalous (i.e. 
inconsistent with our model predictions) polarization at the 
edges of the main jet in its outer parts. We have argued that 
this cannot be a simple superposition of unrelated jet and 
lobe emission. 

(iii) Rapid deceleration is complete before the end of the 
bright region in 3C 296. In contrast, the bright region ends 
within the deceleration zone for B2 0326-1-39, B2 1553-1-24 
and NGC 315 and the end of the bright region roughly co- 
incides with the end of rapid deceleration in 3C31. 

(iv) When deciding where in the jets our hypothesis of in- 
trinsic symmetry first becomes inappropriate, we have suc- 
cessfully used the criterion that the jets must remain straight 
while selecting the region to be analysed. In 3C 296 however 
there are hints of intrinsic asymmetries on a still smaller 
scale because the counter-jet appears to be slightly wider 
than the main jet in the fiaring region. 

Our results are generally consistent with the conclusions 
of a statistical study of a larger sample of FRI sources with 
jets ( Laing et al . 1999), but the velocity range at the bright- 
ening point derived from the distribution of jet/counter- 
jet ratios assuming an isotropic sample has a maximum 
Anax ~ 0.9 and a minimum in the range 0.3 > /?min > 0. 
Such a range would be expected if low edge velocities (as in 
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3C 296) are indeed typical of bridged t win-jet sources, sinc e 
these form the majority of the sample jParma et alj[l99^ . 

8 SUMMARY 

We have made deep, 8.5-GHz images of the nearby FRI ra- 
dio galaxy 3C 296 with the VLA at resolutions ranging from 
0.25 to 5.5 arcsec FWHM, revealing new details of its twin 
jets. In particular we see several thin, discrete brightness en- 
hancements ("arcs") in both jets. These appear to have sys- 
tematically different morphologies in the main and counter- 
jets. A comparison with lower-frequency images from archive 
data shows that the flat-spectrum jets are surrounded by a 
broad sheath of steeper-spectrum diffuse emission, possibly 
formed by a backflow of radiating plasma which has suffered 
significant synchrotron losses as in models of FRII sources. 
The spectral index of the jets initially flattens slightly with 
distance from the nucleus (from a ~ 0.62 to a ~ 0.53) as in 
other FRI jets. We have also imaged the Faraday rotation 
and depolarization over the source. The counter-jet is more 
depolarized than the main jet, but the differences in depo- 
larization and Faraday rotation variance between the two 
lobes are small. The rms fluctuations in rotation measure 
are larger close to the nucleus in both lobes. 

We have shown that many features of the synchrotron 
emission from the inner ±40 arcsec of the jets in 3C 296 can 
be fit accurately and self-consistently on the assumption that 
they are intrinsically symmetrical, axisymmetric, decelerat- 
ing, relativistic flows. The functional forms we use for the 
model are similar to those developed in our earlier work, as 
are most of the derived physical parameters. In this sense, 
a strong family resemblance is emerging between the phys- 
ical properties that we derive for the outflows in these FR I 
sources. 

Geometry The jets can be divided into a flaring region, 
whose radius x is well fitted by the expression x = a2z'^ + 
azz'^ , where z is the distance from the nucleus along the axis, 
and a conical outer region. We infer an angle to the line of 
sight of 58^2 degree. The boundary between the two regions 
is then at 15.8 kpc from the nucleus. 

Velocity structure Where the jets first brighten, their 
on-axis velocity is /3 ~ 0.8, but this value is poorly con- 
strained. Rapid deceleration occurs around 5 kpc to an on- 
axis velocity of /3 = 0.40l:^ []3 which thereafter remains con- 
stant. 3C 296 differs from all of the other objects we have 
studied in having a very low velocity (essentially consistent 
with 0) at the edges of its jets. 

Emissivity The longitudinal profile of rest-frame emissiv- 
ity is modelled as three principal power-law sections e oc 
with indices oi E = 2.5 (0 - 1.8 kpc), 2.8 (2.5 - 8.9 kpc) and 
1.0 (8.9 - 24.5 kpc, the end of the modelled region). The 
first two regions are separated by a short transition zone in 
which the emissivity increases discontinuously by a factor 
of ~7 and then has a power-law slope ~2.1 (1.8 - 2.5 kpc; 
representing the first knots in the main and counter-jets). 
The rest-frame emissivity is centre-brightened. 

Magnetic-field configuration We model the magnetic 
field as disordered, but anisotropic, with toroidal and longi- 
tudinal components, but our results are equally consistent 
with an ordered toroidal component. We rule out a globally- 
ordered helical field. The ratio of longitudinal to toroidal 



field decreases with distance from the nucleus; the toroidal 
component is also more dominant at the edges of the jets. 

Flux freezing The field-component evolution is qualita- 
tively consistent with flux freezing, but much slower than 
expected in the flaring region, where the adiabatic approxi- 
mation also fails to fit the emissivity variation. In the outer 
region, after the jets decelerate, the adiabatic approximation 
is a reasonable fit. 

Arc asymmetry If the arcs are narrow shells of enhanced 
emissivity moving with the underlying flow, then the differ- 
ences in their appearance in the main and counter-jets can 
be understood as the effect of relativistic aberration, as re- 
quired for a fully self-consistent model of the jets. 
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APPENDIX A: FUNCTIONAL FORMS FOR 
VELOCITY, EMISSIVITY AND FIELD 
ORDERING 

In Table \Kl\ we give for reference the functional forms for 
the spatial variations of velocity, emissivity and magnetic- 
field ordering used in this paper. They are identical to those 
used by CLBC except for the elimination of any radial field 
component and a slight change to the form of the longitudi- 
nal emissivity profile. We use the streamline coordinate sys- 
tem (p, s) In the outer region, the streamline index s = S,/^o, 
where ^ is the angle between the (straight) streamline and 
the ax;is. In the flaring region, the distance of a stream- 
line from the jet axis at a distance z from the nucleus is 
x{z, s) = a2{s)z^ + a3{s)z"^, the coefficients being defined by 
continuity of x{s) and its first derivative across the bound- 
ary between the outer and fiaring regions. The coordinate 
p increases monotonically along a streamline and is defined 
by: 



zro sin 
ro sin^o + (cosmos - 1) 
zsin^o + {xo — rosin^o)(l — cosmos) 



sin^o cos ^os 



P < ro 
P^ ro 



On-axis (s = 0), p = 2 in both regions. 
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Table Al. Functional forms of the velocity /3, emissivity e and longitudinal/toroidal magnetic-field ratios k in the streamline coordinate 
system (p, s). Column 4 lists the parameters which may be optimized, for comparison with Table O 



Quantity Functional form 



Range Free parameters 



/3(p, s) /3i exp(-s2 lni;i) 



k{p,s) 



/3icxp( — s In i;i ) — /3q exp( — s lni;o) 
10 



C2 + C3P 



)3o exp(-s^ In i;o) 



+ 



/?! exp( — s lTiv^) — 0Qexp( — s InvQ) 
10 



--Bin 



exp 



Velocity fieW^ 
exp[ci(p - Pvi)] 

exp[c4(pvo - p)] 

Emissivity'' 



, ''"knot 



^2 



exp(— Ineo) 



exp(— In eo) 



p < pvi Distances pvj , Pvq 

Pvi < p ^ Pvo Velocities /3i, /3o 

p > PvQ Fractional edge velocities ni , no 



P ^ P<3i„ Distances pc;,, , pc^^„„^ , Pcout 

Pcin < P ^ Pcknot Indices Ein, -Eknot, i?mid, Eout 

Pcknot < P ^ Pcout Fractional edge emissivities ei, eo 

p > Poout Fractional jump at inner boundary, g 



Longitudinal/toroidal field ratio 



where k^^''' = fcj"''' + (fcg'''^ - fcij"^'^) 



PBn -PBi 



juaxis I c/'i-^ g'^ t.axisN 



P < PBi 

PBi < P < PBo 



P ^ PBo 



Ratios fcf ;,odgc^ ^axis_ ^axis 

Distances pBj , PBq 



" Note that the constants ci , C2 , C3 and C4 are defined by the values of the free parameters and the conditions that the 
velocity and acceleration are continuous at the two boundaries. 

^ The constants d\, and d2 are defined by the condition that the emissivity is continuous except at the inner boundary. 



